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ABSTRACT 


Measurement  of  gas  velocity  in  a  high  enthalpy,  arc-heated  facility, 
has  usually  been  inferentially  derived.  The  ability  to  measure  a  directly 
related  quantity,  the  shift  in  frequency  of  an  emitting  spectral  line, 
enables  the  velocity  to  be  obtained  in  a  noninteracting  manner.  The 
theory  concerning  the  line  shape  and  other  inversion  requirements  is 
discussed.  Three  mathematical  techniques  for  separating  the  lines 
are  discussed  and  examples  of  each  are  given.  Furthermore,  the 
results  of  an  Abel  inversion  are  obtained.  The  test  apparatus  which 
consists  of  the  high  velocity  arc  source,  optics,  spectrometer,  and 
either  a  photographic  plate  or  electro -optical  readout  is  described,  and 
experimental  results  are  presented.  Velocity  measurements  obtained 
by  the  various  experimental  and  mathematical  techniques  are  compared 
as  to  accuracy,  speed,  and  ease  of  data  reduction. 
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SECTION  I 
INTRODUCTION 


During  the  early  1960's  much  emphasis  was  placed  on  the  develop¬ 
ment  of  electric  arc  thrusters  and  magnetohydrodynamic  devices  for 
use  as  space  thrusters.  Simultaneously,  similar  devices  were  developed 
for  use  in  wind  tunnel  applications;  i.  e. ,  reentry,  satellite  and  solar 
wind  simulators.  Most  of  these  devices  produced  environments  hostile 
to  the  use  of  aerodynamic  probe  techniques.  For  example,  high  heat- 
transfer  rates  require  extensive  cooling  of  most  probes.  Aerodynamic 
measurement  of  enthalpy  by  means  of  mass  sampling  techniques  are 
severely  handicapped  by  transitional  flow  effects.  Similar  problems 
also  exist  for  pitot  pressure  probes. 

Although  integrated  measurements  of  introduced  mass  flow,  thrust, 
and  energy  balance  could  be  obtained,  the  average  velocity  could  not. 

This  fundamental  measurement  could  not  be  obtained  directly  from  the 
thrust  equation  because  of  the  large  entrained  mass  flow  values  (Ref.  1). 
For  these  reasons,  a  noninterfering  method  was  sought. 

The  purpose  of  the  study  reported  here  was,  therefore,  to  deter¬ 
mine  the  most  promising  method  by  which  a  routinely  applicable  velocity 
measuring  device  could  be  based.  The  purpose  of  this  report  is  to  sum¬ 
marize  the  approaches  taken  and  to  make  a  critical  comparison  of  differ¬ 
ent  techniques  to  show  the  direction  in  which  a  final  development  of  such 
a  device  will  have  to  go. 


SECTION  II 
THEORY 


Measurement  of  the  Doppler  shift  of  a  spectral  emission  line  re¬ 
quires  a  knowledge  of  the  line  profile.  The  following  sections  discuss 
these  effects.  The  modification  of  spectral  lines  attributable  to  instru¬ 
mental  broadening,  although  admittedly  important,  has  not  been  included. 
Further  work,  which  includes  inversion  of  the  observed  lines,  should 
take  these  effects  into  consideration. 
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2.1  DOPPLER  SHIFT  AND  LINE  SHAPE 

The  apparent  change  in  wavelength  of  radiation  caused  by  a  velocity 
difference  between  the  radiation  source  and  observer  is  given  by 

-  X'  =  A0  V/c 

Where  V  is  taken  to  be  positive  when  the  source  is  moving  toward  the 
observer,  A0  is  the  unshifted  wavelength,  and  A'  is  the  shifted  wave¬ 
length.  This  shift  can  be  experimentally  obtained  by  observing  the  radi¬ 
ation  from  two  different  optical  paths  through  the  luminous  flow  as  in 
Fig.  1.  For  a  uniform  flow  velocity  V,  the  shift  in  wavelength  between 
the  two  lines  is 

,  /  v 

Ad  =  Ai  -  A2  =  A0  —  (cos  ©i  +  cos  62 ) 

Again  the  primed  symbols  indicate  the  wavelength  of  the  shifted 
lines.  (For  clarity,  incremental  symbols  have  been  dropped  from  Ad.) 
Although  any  combination  of  angles  could  be  chosen,  they  were  made  as 
nearly  equal  as  possible  to  minimize  effects  caused  by  axially  dependent 
variation  in  velocity  or  intensity.  Of  course,  the  angle  should  be  mini¬ 
mized  to  maximize  the  measurable  shift  and  thereby  improve  accuracy. 

Spectra  were  recorded  either  directly  on  photographic  plates  or  on 
an  oscilloscope  after  first  being  detected  by  a  photomultiplier  tube  and 
scanning  device.  In  either  case,  it  was  possible  to  record  the  radiation 
from  both  lines  in  two  ways:  either  simultaneously,  which  resulted  in 
an  unresolved  composite  line,  or  at  two  different  times  which  lead  to 
distinct  lines  (Section  2.  3). 

When  the  lines  are  distinct,  an  average  velocity  can  be  obtained  by 
measuring  the  relative  peak  displacement  with  little  regard  for  the  line 
shape.  However,  for  unresolved  lines,  the  shape  of  the  composite  line, 
when  compared  to  that  of  a  single  distinct  line,  contains  the  Doppler 
shift  information. 

The  line  shape  for  dipole  radiation  of  a  single  spectral  line  is  given 
by  the  classically  derivable  Lorentzian  line  profile  when  written  in 
terms  of  frequency 
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where 

,  87r2e2ro2  ,  .  _  f  , 

y  =  — gmc3 —  (cgs  units);  IQ  =  J  I(y)di/ 

The  natural  line  width  is  characteristically  small:  1.  16  x  10“^ A.  How 
ever,  in  the  case  of  thermally  energetic  gases  the  thermal  motion  can 
broaden  line  radiation  into  the  Doppler  profile  which  is  Gaussian.  The 
relative  intensity  as  a  function  of  frequency  is 

Kv)  =  IDe 

The  full  width  at  half  of  maximum  intensity  is 

1.  67  —  1.67^2  IE T 

1 'oyj  m  C  yj  m 

This  relation  is  frequently  used  to  determine  the  heavy  particle  temper¬ 
ature. 

Figure  2  shows  a  Gaussian  curve  fit  to  the  4266.  53 A  Ar  II  line 
obtained  by  a  photographic  method  and  a  photoelectrical  method.  Each 
has  been  arbitrarily  normalized.  The  difference  in  their  half-widths  is 
attributable  to  instrumental  broadening  effects.  In  the  photographic 
method  the  line  was  reduced  from  the  plate  by  means  of  a  microphotom¬ 
eter  and  strip- chart  readout,  whereas  the  photoelectric  line  was  taken 
directly  from  an  oscilloscope  trace.  The  selection  of  an  accurate  base 
line  is  critical  for  certain  methods  as  will  be  shown  later.  Actually,  a 
Gaussian  curve  should  not  in  general  be  a  good  fit  as  the  line  is  a  sum 
over  the  entire  optical  path.  For  example  a  skewed  line  can  result  even 
though  the  contributing  effects  are  gaussian  (Fig.  3).  Of  course,  the 
flow  properties  are  continuous  rather  than  discrete,  and  it  is  the  function 
of  an  inversion  technique  to  include  these  effects. 


2.2  ABEL  INVERSION 

The  line  intensity  I(x)  is  related  to  the  radial  emission  coefficient 
I(r)  in  the  usual  manner  by  the  relation  (Refs.  3,  4,  5,  and  6) 

t/  \  I(r)  rdr 

I(x)  =  2Jx-p^j1/2 
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Fig.  3  Distortion  of  Composite  Line 
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The  inversion  is  given  by 


I(r) 


.1  I  I  (x)  dx 
?r'r-(x2  -r2)1/2 


v.  _  dl 

where  I  =  3— 
dx 


However,  in  the  case  of  finite  station  widths  (Fig.  4)  an  incremental 
approach  must  be  taken.  The  plume  is  considered  to  be  composed  of  a 
series  of  concentric  rings  of  unit  spacing  and  unit  depth.  This  assump¬ 
tion  is  valid  even  though  the  angle  of  view  is  not  normal  to  the  flow  axis. 
The  viewing  angle  increases  the  observed  intensities,  I(x)'s,  by  the 
multiplier  esc  9  for  the  cylindrically  symmetric  case.  The  intensity  ob¬ 
served  at  station  i  is  assumed  to  be  a  sum  over  all  the  rings  contained 
in  the  slice,  weighted  by  elemental  areas.  That  is 


J  =n 

h  =2  hJi 

j  =  l 


In  matrix  form 


where 


A 

I  =  A  I 


V 

h*  — 

h 

I  = 

12 

■ 

• 

A 

...  I  = 

A 

12 

In 

A 

A  =  jAi|A2|...|Anl 


The  A's  which  form  the  triangular  matrix  A  are  found  in  the 
Appendix. 


A  1 

To  determine  the  radial  distribution  one  simply  forms  I  =  A  I. 
As  A  is  triangular,  Aa  Gauss  reduction  can  be  performed  yielding 
rather  directly  the  I's, 

Since  only  relative  intensity  was  required,  units  were  not  rigidly 
adhered  to.  Any  scaling  was  performed  after  inversion  and  was 
merely  as  a  matter  of  convenience. 
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Fig.  4  Abel  Inversion 
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Input  data  were  composed  of  n  profiles  of  I  versus  X  with  XQ  corre¬ 
sponding  to  Xd  =  0  obtained.  This  provided  a  three-dimensional  array 
as  shown  in  Fig.  5.  The  inversion  program  enters  at  a  constant  X  and 
inverts  the  data.  The  number  of  such  inversions  depends  on  the  number 
of  data  points  making  up  each  input  profile;  the  three-dimensional  array 
of  I,  r,  X  is  formed  (Fig.  6).  These  are  the  inverted  data,  and  their 
treatment  is  described  below. 


2.3  DOPPLER  SHIFT  MEASUREMENTS 

The  measurement  is  either  in  the  form  of  a  composite  unresolved 
line  plus  a  single  line  or  two  distinct  lines  recorded  at  different  times. 
The  exact  form  determines  the  means  of  data  reduction. 

When  the  lines  are  distinct,  measurements  can  be  made  in  one  of 
two  ways: 

1.  If  the  two  differently  shifted  lines  can  be  independ¬ 
ently  superimposed,  such  as  they  could  in  a  photo¬ 
electric  recording,  the  measurement  of  the  shift 
can  be  made  directly  by  measuring  the  peak  dis¬ 
placements. 

2.  When  the  two  lines  were  on  different  spectra,  as 
with  the  photographic  plate,  a  reference  system  of 
lines  from  a  standard  source  was  superimposed. 

Measurements  were  then  taken  with  respect  to  a 
nearby  standard  reference  line.  In  both  cases  two 
reference  lines  are  required  to  obtain  the  system 
dispersion. 

When  noise  is  a  factor,  whether  caused  by  electronic  or  emulsion 
grain  effects,  several  estimates  of  the  location  of  the  center  can  be 
made.  The  center  of  an  intensity  trace  can  be  estimated  several  times 
by  obtaining  the  midpoint  of  the  full  width  at  various  fractions  of  the 
maximum  line  intensity.  The  center  of  a  line  viewed  directly  from  the 
photographic  plate,  as  with  a  contour  projector,  can  be  estimated  by 
taking  several  measurements  along  the  length  of  the  line. 

The  Doppler  shift  can  be  obtained  from  the  characteristics  of  a 
composite  and  single  line  shape.  Three  such  techniques  were  developed, 
each  requiring  a  different  form  of  the  data.  As  the  algebraic  manipu¬ 
lations  are  straightforward,  only  the  general  form  and  the  results  for 
a  Gaussian  line  shape  are  shown. 


9 


AEDC-TR-71-175 


Most  line  profiles  are  given  in  the  form  I  =  IQ  f(X  -  a)  where  (a)  is 
assumed  to  be  the  center  of  the  line  and  corresponds  to  that  wavelength 
where  I  =  ID.  The  composite  profile  consists  of  the  sum  of  two  equal 
intensity  profiles  which  are  shifted  from  each  other  by  an  amount  Xq, 
the  Doppler  shift.  If  all  measurements  are  made  relative  to  the  center 
of  this  line,  then  the  composite  profile  is  given  by 


Ic<*>  *  Iq 


The  first  method  uses  a  parameter  a  which  is  the  ratio  of  the 
composite  line  intensity  evaluated  at  its  center  to  that  of  a  single  line 
at  its  center;  i.  e. , 


The  Doppler  shift  for  a  Gaussian  distribution,  when  normalized  to  the 
single  line's  half -width  at  one  half  maximum  intensity,  6S,  is 


Figure  7  indicates  that  for  large  separations  the  method  becomes  very 
sensitive  to  small  changes  in  intensity  ratios.  The  main  shortcoming  of 
this  method  is  that  great  care  must  be  taken  when  exposures  are  made. 

A  very  slight  imbalance  attributable  to  exposure  time,  for  example,  can 
cause  considerable  error.  The  instrument  as  utilized  could  not  readily 
achieve  this  criterion  as  all  shutter  operations  were  manually  controlled. 


The  second  method  utilizes  the  parameter  3  which  is  the  ratio  of 
the  intensity  of  the  composite  line,  evaluated  at  ±  6S  from  its  center, 
to  the  composite's  intensity  at  its  center;  i.  e. , 

B±  ntw  _f(**s  +  irH(±{s-ir) 

=  Mo>  '  ¥)+<-¥) 

Again  the  Doppler  shift  for  a  Gaussian  distribution  is 

Xp  =  tn  (23  ±^432  -  l) 

6s  Pn  2 

This  is  shown  plotted  in  Fig.  8.  Since  measurements  are  taken 
from  a  single  line  profile  composite  and  the  invariant  half-width, 
exposure  matching  is  not  a  problem. 
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The  final  method  utilizes  the  ratio  of  the  half-widths  of  a  single  and 
a  composite  line  as  measured  at  one  half  the  center  intensity.  This 
ratio  is 


With  no  shift  this  ratio  is  unity  and  increases  with  increasing  A.p).  This 
method  requires  only  the  measurement  of  the  half -widths  to  obtain  the 
shift  and  perhaps  is  the  most  straightforward  approach. 


The  general  implicit  relation  is 

IC(±6C)  =  [M°>]/2. 


For  the  Gaussian  line  shape  (Fig.  9) 

(  72-l  /  2<t2-1)  \ 

_  en\2  ±V2 _ -1/ 

6s  Tfn  2 


SECTION  III 
EQUIPMENT 


The  apparatus  consisted  of  the  high  velocity  arc  source,  inter¬ 
mediate  optics,  a  dispersion  element  in  the  form  of  a  spectrometer, 
and  either  a  photographic  plate  or  a  scanner  and  a  photomultiplier  tube. 
Photographic  data  reduction  required  a  photodensitometer  or  a  traveling 
microscope.  When  data  were  to  be  taken  electronically,  an  oscillo¬ 
scope  triggered  by  the  scanner  was  used  to  make  photographic  record¬ 
ings  of  the  data. 


3.1  HIGH  VELOCITY  ARC  SOURCE 

The  tests  were  conducted  using  a  self -magnetic  arc  accelerator. 

The  electrodes  for  the  d-c  arc  were  configured  to  utilize  the  self- 
induced  magnetic  field  of  the  discharge  to  accelerate  the  ionized  argon 
working  gas  (Fig.  10).  The  current  flows  from  the  anode  to  the  cathode. 
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heating,  ionizing  and  accelerating  the  gas  as  it  expands.  The  plume  is 
diffuse  with  no  evidence  of  "spoke"  formation.  At  the  point  of  observa¬ 
tion,  just  downstream  of  the  anode  ring,  the  plume  was  axisymmetric 
and  about  6  in.  in  diameter. 

The  device  was  contained  in  a  4-  by  4-  by  8-ft  vacuum  chamber 
(Fig.  11).  Optical  access  was  provided  by  quartz  windows  in  the  test 
cell. 


A  combination  of  pumps  was  used  to  obtain  typical  test  cell  pres¬ 
sures  of  100  A*Hg. 

The  Eire  power  supply  was  current  controlled  and  capable  of 
supplying  current  up  to  6000  amp.  All  data  were  taken  at  the  following 
conditions : 

Current  1500  amp 

Voltage  53.  5  volts 

Working  Gas  Argon 

Mass  Flow  Rate  0.075  gm/sec 

Although  the  arc  parameters  remained  quite  constant,  the  light  intensity 
from  a  single  spectral  line  decreased  in  time  (Fig.  12).  This  phenom¬ 
enon  paralleled  the  heating  of  engine  components  and  could  be  caused  by 
quenching  of  metastable  states  by  impurities.  .  No  effect  was  observed 
upon  the  velocity  measured  at  different  times  into  a  run,  utilizing  the 
peak-peak  shift  method.  Since  data  acquisition  times  were  not  coinci¬ 
dent,  a  systematic  error  would  result  when  other  data  reduction  tech¬ 
niques  were  used,  because  intensity  matching  was  important.  To  reduce 
this  problem,  data  were  taken  about  60  sec  after  arc  startup. 


3.2  SPECTROMETER 

A  Jarrell  Ash  1.  5-m  grating  spectrometer  Model  75-000  was  used. 

To  avoid  vibration,  the  spectrometer  was  isolated  from  the  test  cell  and 
pumps  by  a  separate  stand  with  foam  rubber  pads.  Although  small  dis¬ 
placements  could  exist  between  the  arc  and  the  spectrograph  because  of 
the  mounting  arrangement,  this  effect  was  considered  to  be  negligible. 

All  optical  components  were  mounted  either  directly  to  the  spectrometer's 
optical  bench  or  to  a  0.  5-in.  -thick  aluminum  plate  secured  to  it. 

Slit  width  capabilities  are  from  4  to  400 p,  where  the  lower  limit  is 
controlled  by  an  automatic  stop  to  protect  the  edges.  The  practical  lower  - 
limit,  because  of  diffraction,  is  about  10 M.  To  shorten  required  run 
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times  the  slit  widths  used  were  increased  from  17  to  20  This  had 
little  effect  on  tjie  line  widths;  e.  g. ,  a  3-M0slit  adjustment  would  in¬ 
crease  a  0.  20-A-wide  Ar  II  line  to  0.  215  A.  When  scanning  onto  the 
phototube,  a  compromise  was  made  between  resolution  and  the  signal  - 
to-noise  ratio.  Slit  height  has  only  second-order  effects  on  the  reso¬ 
lution.  A  3 -mm  height  was  used  throughout  the  tests. 

The  spectrometer  also  had  to  be  adjusted  for  maximum  resolution. 
This  was  achieved  by  adjustments  at  both  the  entrance  slit  and  the 
focal  plane  at  the  exit  of  the  spectrometer  by  use  of  the  iron  arc  source. 
Adjustments  were  first  made  for  the  photographic  plate  by  a  series  of 
exposures  at  different  entrance  slit  widths  and  plate  angles.  When 
this  was  finished,  the  photomultiplier  assembly  and  the  exit  slit  were 
inserted  and  optimized  by  adjusting  the  assembly  and  rotating  the 
entrance  slit.  In  this  manner  both  the  photographic  and  photoelectric 
assemblies  were  optimized  and  thereafter  could  be  interchanged  at  will. 


3.3  STANDARD  SOURCE 

A  standard  source  for  wavelength  calibration  was  a  Jarrell-Ash, 
Fe-Ne  hollow  cathode  arc  which  had  a  maximum  current  of  20  ma.  The 
d-c  power  supply  was  regulated  and  had  a  ripple  factor  of  0.  05  percent. 
The  source  was  mounted  on  the  spectrometer  where  it  could  be  focused 
on  the  entrance  slit  by  means  of  a  lens  and  a  flip  mirror. 


3.4  INTERMEDIATE  OPTICS 

The  purpose  of  the  intermediate  optics  was  to  collect  radiation  from 
the  plume  along  two  predetermined  paths,  combine  them  in  such  a 
manner  that  they  were  collinear,  and  focus  this  radiation  on  the  entrance 
slit  of  the  spectrometer.  Furthermore,  symmetry  conditions  were  to 
be  maintained  while  the  maximum  possible  velocity  component  was  to  be 
viewed. 

The  system  functions  as  follows  (Fig.  1).  Light  from  the  two  paths 
passes  through  the  windows  W^  and  W2  and  is  directed  to  the  beam 
splitter,  BS,  by  mirrors  Mi  and  M2.  Here  they  are  combined  and 
directed  through  the  aperture  Ai  and  the  lens  L  by  the  mirror  M5.  The 
lens  focuses  the  arc  centerline  image  on  the  spectrometer  entrance 
slit,  Sj.  Auxiliary  movable  mirrors  M3  and  M4  direct  the  iron  arc  and 
the  helium-neon  alignment  laser,  respectively. 
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All  mirrors  were  front -surfaced  and  aluminum  coated.  The  beam 
splitter  was  aluminized  quartz  with  with  a  nominal  40-percent  trans¬ 
mission  coefficient. 

Aperture  Aj  serves  primarily  as  a  field  stop  to  maintain  an 
f/number  of  12.  This  prevents  overfilling  the  spectrometer  focusing 
mirrors.  The  aperture  was  located  as  dose  to  the  lens  as  possible  so 
that  the  relative  weighting  of  the  contributions  along  the  optical  path  in 
the  plume  was  unchanged  by  changes  in  aperture  setting.  When  only  one 
set  of  measurements  was  to  be  made,  the  system  focused  on  the  arc. 
centerline  rather  than  at  infinity.  This  weighed  all  portions  of  the  plume 
equally  as  the  optical  and  plume  geometries  were  nearly  congruent. 

3.4.1  Alignment 

Alignment  of  the  intermediate  optics  was  achieved  through  the  use 
of  a  laser  and  a  target  system.  Mirror  M4  was  placed  in  front  of  the 
entrance  slit  and  a  removable,  prealigned  set  of  targets  was  placed  on 
the  optics  bench.  Minor  adjustments  of  the  laser  and  mirror  provided 
a  reference  beam  which  was  on  the  spectrometer  optical  axis.  From 
this  point  adjustments  were  made  working  backward  to  the  two  mirrors 
Mi  and  M2.  Two  very  small  pieces  of  reflective  tape  on  the  far  test 
cell  wall  were  used  as  targets  to  achieve  this  final  alignment.  Consider¬ 
able  care  was  taken  tc>  eliminate  or  minimize  all  other  reflected  light 
through  blacking  of  the  test  cell  surfaces.  The  results  of  improper  pre¬ 
cautions,  where  a  corner  reflector  was  inadvertently  produced,  are 
discussed  later. 

Rather  than  calculate  the  focal  points  for  the  lens,  another  lens  was 
placed  in  front  of  the  laser  to  focus  the  beam  to  a  point.  The  distance 
from  this  point  to  mirror  M4  was  adjusted  to  coincide  with  that  from  the 
mirror  to  the  entrance  slit.  Si.  Then  by  turning  either  Mi  or  M2  so  that 
the  beam  fell  outside  the  test  cell,  Li  was  moved  so  that  the  image  of 
the  point  was  a  minimum  at  the  premeasured  distance.  Corrections 
were  not  made  for  a  very  slight  distortion  caused  by  the  quartz  windows 
because  the  point  of  the  beam  crossover  was  predetermined  geometri¬ 
cally.  . 

3.4.2  Beam  Balance 

The  actual  balancing  of  the  two  beams  was  achieved  by  either  vary¬ 
ing  an  aperture  in  one  of  the  arms  or  by  adding  neutral  density  filters. 

The  former  method  causes  a  change  in  the  relative  distribution  over  the 
elements  of  the  optical  path.  The  latter  is  limited  because  of  the  cutoff  of 
the  glass  substrates  at  about  3900  A.  By  using  quartz,  this  limit  can  be 
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moved  down  to  <  2000  A.  This  transmission  characteristic  causes  a 
variation  in  intensity  over  the  entire  recorded  spectrum.  The  net  result 
is  an  intensity  match  only  at  the  monitored  wavelength.  This  nullifies 
to  some  extent  the  advantage  of  the  photographic  methods,  that  of  simul¬ 
taneously  obtaining  a  large  quantity  of  information.  This  problem  can 
be  eliminated  through  the  choice  of  a  beam  splitter  and  the  elimination 
of  other  nonsymmetric  frequency-dependent  attenuators  such  as  extra¬ 
neous  coating  on  optics.  For  example,  the  quartz  viewing  ports  became 
heavily  coated  with  various  copper  compounds  after  several  hours  of 
arc  running  which  cause  preferential  attenuation. 


3.5  MICROPHOTOMETER 

The  recording  microphotometer  is  designed  for  precision  measure¬ 
ment  of  the  blackening  of  photographic  emulsion.  It  magnifies  by  15X 
the  image  of  the  plate  utilizing  a  high  illumination  source.  As  a  photo¬ 
graphic  plate  is  scanned  the  images  of  the  spectral  lines  move  across  a 
slit;  the  intensity  of  the  images  is  proportional  to  the  light  transmitted 
by  the  plate.  Thus  the  current  output  of  its  photomultiplier  is  propor¬ 
tional  to  the  transmission  coefficient  of  the  line.  The  current  is  ampli¬ 
fied  and  used  to  drive  a  strip-chart  recorder.  Thg  microphotometer  is 
capable  of  speeds  as  low  as  0.  5  mm/min  or  0.  25  A/min. 

The  carriage  motion  was  insufficiently  reproducible  over  long  dis¬ 
tances  on  the  plate  because  of  motor  and  gear  slack.  Measurements 
taken  from  the  microphotometer  yielded  essentially  the  same  results 
as  those  of  the  contour  projector  (3.  6  below);  however,  many  more 
passes  were  required  to  approach  the  same  level  of  confidence. 


3.6  CONTOUR  PROJECTOR 

A  contour  projector  was  also  used  to  read  displacements  from 
photographic  plates.  Measurements  to  0.  0001  in.  can  be  made  (which 
corresponds  to  0.  0127  A).  The  projector  allowed  manual  reading  of 
carriage  displacements  as  viewed  from  a  projection  screen.  Magnifi¬ 
cations  from  10X  to  100X  were  possible,  and  the  lowest  possible 
magnification  compatible  with  measurement  accuracy  was  used.  The 
contour  projector  was  most  satisfactory,  especially  since  it  permitted 
human  judgment  which  was  not  inherent  in  the  microphotometer. 
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SECTION  IV 

PHOTOELECTRIC  SCANNING 


The  equipment  required  for  the  method  is  shown  in  the  block 
diagram  (Fig.  13).  Light  from  the  spectrometer  was  scanned  across 
a  10-p  exit  slit  onto  the  photomultiplier  tube.  This  was  in  turn  ampli¬ 
fied  and  filtered  before  entering  the  oscilloscope.  Display  synchro¬ 
nization  was  obtained  through  a  trigger  circuit.  Certain  optional  com¬ 
ponents  are  indicated  but  were  not  utilized. 


4.1  REFRACTOR  PLATE 

The  scanning  was  accomplished  by  means  of  a  rotating  refractor 
plate  located  in  front  of  the  exit  slit.  A  small  d-c  motor  was  used  to 
turn  the  refractor  plate  at  speeds  as  low  as  2  rpm  (Fig.  14).  Also 
shown  in  the  photo  is  the  microswitch  and  cam  which,  through  a  filter 
circuit,  triggered  the  oscilloscope.  A  leaf  spring  was  added  (not  shown) 
which  applied  pressure  against  the  cam  collar,  reducing  backlash  in 
the  gear  train.  The  amount  of  displacement  of  the  light  rays  for  a 
given  angular  displacement  was  assumed  to  be  linear  for  near-normal 
incidence. 


4.2  ELECTRONICS 

Figure  15  shows  the  electronics  schematic.  The  photomultiplier 
tube  (10-mm-diam  photocathode)  was  supplied  by  a  2000-V  dc  power 
supply.  The  signal  from  the  photomultiplier  tube  passes  through  an 
input  filter  consisting  of  a  1-Mfi  resistor  in  series  and  a  0. 002  Mf 
capacitor  in  parallel  for  noise  elimination.  The  filtered  signal  then 
enters  a  variable  bandpass  filter  preamplifier  plug-in  unit  (response 
from  0  to  100  Hz)  where  further  noise  elimination  takes  place.  The 
signal  is  then  displayed  on  the  oscilloscope  where  a  conventional  scope 
camera  was  used. 

The  triggering  unit  removes  voltage  pulses  attributable  to  switch 
chatter,  eliminating  false  triggering.  There  is  no  negative  pulse  upon 
release  of  the  switch  because  of  the  delay  constant  of  the  low  pass 
section  (0.  06  sec).  The  high  pass  section  does  not  respond  to  this 


23 


Fig.  13  Photoelectric  Scanning 
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delay.  Variations  in  C2  and  R2  can  be  made  to  match  input  impedence 
if  it  is  not  50  fi.  With  a  50  ^  load,  R2  was  not  required. 


SECTION  V 
DATA  ACQUISTION 


After  the  intermediate  optics  and  spectrometer  were  aligned,  a 
prerun  check  was  made.  With  the  photoelectric  equipment  in  place  and 
with  the  scanner  running  at  a  nominal  2  rpm,  the  two  Fe  lines,  4271. 16 
and  4271.  76  A,  were  observed  for  a  calibration,  intensity,  and  focus 
check.  Also,  if  the  signal -to-noise  ratio  was  poor,  corrective  steps 
were  taken.  The  arc  was  then  started  and  allowed  to  stabilize.  During 
this  time  the  Ar  II  line  to  be  used  was  located  and  appropriately 
centered  on  the  oscilloscope.  This  was  achieved  with  the  aid  of  the 
grating  angle  adjustment. 


5.1  INTENSITY  MONITOR 

Variations  in  the  intensity  level  of  the  two  lines  over  an  extended 
series  of  shots  resulted  in  data  which  were  considered  unsuitable  for 
use  with  the  Abel  inversion  program.  Therefore,  a  monitoring  system 
was  arranged  and  tested.  It  consisted  of  a  0.  25-m  spectrometer  and 
associated  optics.  It  viewed  the  same  axial  station  as  did  beams  1 
and  2.  Its  output  was  placed  on  a  dual -be  am  oscilloscope  so  that  a  moni 
tor  level  and  signal  trace  were  recorded  simultaneously.  After  the  data 
were  analyzed,  the  level  was  used  to  normalize  the  line  profiles.  The 
results  of  this  study  indicated  a  variation  of  less  than  7  percent  while 
data  were  being  obtained.  Corrections  based  on  this  monitor  system 
did  not  improve  the  data.  However,  it  cannot  be  said  that  such  a  tech-  • 
nique,  especially  if  used  as  an  automatic  gain  control,  would  not  con¬ 
tribute  substantially  to  the  accuracy  of  the  Abel  data. 

For  general  application,  however,  a  line  intensity-time  profile  is 
desirable,  especially  if  the  time  of  measurement  is  on  the  order  of  the 
time'  required  for  a  nominal  5  percent  level  change  to  occur. 


5.2  PHOTOGRAPHIC  DATA  ACQUISTION 

The  type  of  photographic  plate  used  to  record  the  data  was  a  com¬ 
promise  between  resolution,  rim  time,  and  wavelength  response. 


27 


AEDC-TR-71-175 


The  103 F  plate  had  good  sensitivity  from  the  ultraviolet  up  to  about 
5700  A  with  medium  grain  size  and  moedium  contrast.  The  slower  speed 
SAI  plates  had  a  cutoff  of  about  4300  A,  with  fine  grain,  medium  con¬ 
trast,  and  higher  resolution.  Because  of  the  optics  cutoff  wavelength 
(3900  to  2000  A),  below  which  the  light  was  readily  attenuated,  the 
range  ofoavailable  wavelength  was  limited  to  not  much  more  than  4200 
to  4300  A  for  glass  or  2000  to  4300  A  for  quartz  with  the  SAI  plate. 

The  corresponding  ranges  for  the  103F  plate  were  greater— 3900 
to  5700  A  for  glass  and  2000  to  5700  A  for  quartz  optics. 

Exposure  times  for  the  103  F  using  a  20 -M  entrance  slit  were  0.  5 
to  5  sec  for  the  argon  lines  and  15  to  40  sec  for  the  iron  arc  source 
medium  to  strong  lines. 

Background  fogging  was  found  to  be  high  if  the  plates  were  not 
warmed  from  cold  storage  temperatures  before  use. 

Exposures  taken  over  extended  periods  of  time  required  dimen¬ 
sional  stability  of  the  plate  and  camera  assembly,  which  required  an 
elaborate  method  of  spectrometer  isolation. 

Data  were  taken  in  the  manner  indicated  in  Table  I.  Each  run  con¬ 
sists  of  various  combinations  of  the  radiation  taken  from  either  the 
upper  stream,  downstream  arc  source  or  the  iron  arc  standard.  Rack 
positions  correspond  to  different  run  numbers  and  to  different  exposure 
bands  on  the  plate. 


SECTION  VI 
DATA  REDUCTION 


Data  reduction  methods  were  largely  a  function  of  the  method  of 
acquisition.  The  photoelectric  data  could  be  directly  reduced  by  means 
of  the  peak  displacement  method.  Also,  if  a  composite  were  taken,  the 

а,  j3,  and  y  (or  line  shape)  methods  could  be  applied.  For  the  photo¬ 
graphic  plate  method,  the  photographic  plate  first  had  to  be  developed, 
then  the  data  were  either  reduced  to  an  intensity  trace  by  means  of  a 
photodensitometer,  or  relative  displacements  were  measured  by  the 
contour  projector.  In  the  case  of  the  photodensitometer  method,  the  a, 
3,  and  y  methods  were  applicable. 

б. 1  PHOTOELECTRIC  REDUCTION 

Electronic  data,  when  taken  as  a  photograph  of  the  oscilloscope 
trace,  were  immediately  ready  for  data  reduction.  Peak-to-peak 
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TABLE  I 

PHOTOGRAPHIC  PLATE  IDENTIFICATION 


Run 

Rack 

Time, 

Beam 

Direction,  Source 

Slit  Height, 

Slit  Width, 

No. 

No. 

Sec 

mm 

1 

300 

5 

1,  Ar  II 

* 

2 

0 

2 

350 

5 

2,  Ar  II 

3 

400 

2.5 

1  and  2,  Ar  II 

A 

500 

5 

1,  Ar  II 

4 

15 

Fe  Arc 

5 

550 

5 

15 

2,  Ar  II 

Fe  Arc 

6 

600 

2.5 

1  and  2,  Ar  II 

15 

Fe  Arc 

*Beam  direction  (1  or  2)  is  described  by  the  angles  0^  and  62 
in  Fig.  1. 


Plate  Type:  103F  Emulsion 
Grating  Angle:  14.  00  deg 
Date:  1/1/70 
Data  Point  Number:  200 
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measurements  were  made  directly  from  the  traces  as  shown,  in  Fig.  16. 
When  noise  was  a  factor,  several  estimates  of  peak  X  were  made  by  . 
graphically  deriving  the  half-width  of  the  envelope  at  selected  fractional 
values  of  IQ.  Unless  one  of  the  lines  was  badly  skewed  by  reflections  of 
oppositely  shifted  radiation  within  the  test  cell  (corner  reflection)  dis¬ 
placements  obtained  in  this  manner  agreed  with  direct  peak  measure¬ 
ments.  Calibration  was  obtained  by  simultaneously  scanning  the  iron 
lines,  which  are  separated  by  0.  605  A  (Fig.  17).  This  is  not  the  only 
possible  method  for  obtaining  the  dispersion.  An  iron  line  and  even  a 
shifted  Ar  II  could  have  been  used  if 


*Fe  "  *Ar  I 

For  most  lines  this  criterion  is  easily  met.  'J’he  average  shift  was 
0.  0608  ±  0. 0063  A  with  a  resolution  of  0. 103  A.  For  the  viewing  angles 
of  64. 1  deg,  this  corresponds  to  a  velocity  of  4340  ±  474  km/sec. 


As  a  check,  the  a,  j3,  and  7  methods  were  used  to  determine  their 
reproducibility  and  limitations  because  of  the  approximations  used. 

For  all  these  methods  a  base  line  was  drawn  which  approximately 
averaged  out  the  noise  level.  As  a  reasonable  curve  fit  was  obtained, 
the  Doppler  shift  values  based  on  the  Gaussian  curve  were  used.  Mis¬ 
matches  in  intensity  were  averaged  whenever  two  measurements  of  the 
same  quantity  could  be  obtained  for  a  single  variable;  e.  g. ,  when  the  7 
method  is  used,  two  measurements  of  the  two  single  line  half -widths 
could  be  averaged.  The  results  are  tabulated  in  Table  II. 


6.2  PHOTOGRAPHIC  REDUCTION 

The  microdensitometer  comparator  with  a  strip-chart  recorder  out¬ 
put  was  used  originally  to  obtain  a  trace  of  the  line  profiles.  Although 
the  accuracy  of  the  instrument  was  adequate  for  intensity  studies,  a 
wavelength  error  of  ±0.  025  A/in.  was  not  sufficiently  small  to  measure 
differences  from  0.  02  to  0.  07  A.  Measurements  were  made  using  this 
device,  but  a  large  quantity  of  redundant  densitometer  tracings  were 
required  for  each  line  and  are  not  included  here. 

Data  tabulated  in  Table  III  were  taken  with  the  contour  projector. 

AH  measurements  were  taken  relative  to  fixed  Fe  lines,  but  were 
numerically  compared  to  the  composite  Ar  II  line.  This  artifice 
served  no  purpose  other  than  to  indicate  mismatches  in  intensity.  The 
net  shift  is,  therefore,  the  algebraic  stun  of  these  differences. 
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Fig.  16  Shifted  Argon  Lines 


Fe  -1271.764  K 


4271.159  X 


Fig.  17  Iron  Calibration  Lines 
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TABLE  II 

TEST  RESULTS  OF  PHOTOELECTRIC  TECHNIQUE 


Wavelength  4266  A 


Reduction  Method 

Peak-to-Peak 

a 

0 

7  | 

A\d>  a 

AXd,  A 

AAd 

AAd 

0.0579 

0.0958 

0.0614 

0.0211 

0.0579 

0.0896 

0.0806 

0.0719 

0.0527 

0.0780 

0.0850 

0.0795 

0.0754 

0.0186 

0.0680 

0.0675 

0.0560 

0.0654 

0.0720 

0.0763 

0.0615 

0.0834 

0.0673 

0.0930 

0.0560 

0.0630 

0.0763 

0.0527 

0.0707 

0.0682 

0.0670 

0.0544 

0.0598 

0.0688 

0.0755 

0.0759 

0.0688 

0.0783 

0.0773 

0.0600 

0.0695 

0.  0744 

0.0670 

0.0581 

0.0554 

0.0583 

0.0565 

0.0565 

0.0666 

0.0765 

Average 

0.0608  ±0.0066 

0.0718 

0.0693 

0.0695 

V 

4340  ±  474 

5140  ±  2010 

4960  ±  400 

4970  ±  1180 
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TABLE  III 

TEST  RESULTS  OF  PHOTOGRAPHIC  TECHNIQUE 


Wavelength  (XQ),  A 

(AA.d)Mo  x 

4426 

0.020 

4.53 

4228 

.  0. Ill 

25.9 

4228 

0.046 

10.  7 

4228 

0.071 

16.5 

4013 

0.039 

9.4 

4013 

0.023 

5.7 

4013 

0.076 

18.  75 

3928 

0.  113 

27.  8 

3850 

0.050 

13.0 

3850 

0.091 

23.  6 

3850 

0.067 

17.4 

3765 

0.040 

10.  65 

3765 

0.080 

21.3 

3729 

0.  087 

23.2 

3729 

0.036 

9.  65 

3588 

0.067 

17.25 

3588 

0.045 

12.  55 

3576 

0.073 

20.4 

3576 

0.017 

12.8 

Average:  15.  84  ±  6.  72 

Velocity:  5400  ±2310  m/sec 

la  on  mean  530  m/sec 
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6.3  ABEL  INVERSION 

The  results  of  the  Abel  inversion  are  shown  in  Figs.  18  and  19. 
Note  that  a  strong  peaking  occurs  at  ring  5.  This  is  believed  to  be 
caused  by  the  presence  of  barrel  shocks  within  the  plume  itself.  As 
can  be  readily  observed,  actual  peak  location  is  in  question.  Several 
methods  were  used  to  determine  the  peak:  (1)  determination  of  the 
centroid,  (2)  peak  displacement,  and  (3)  a  series  of  midpoint  values. 
None  gave  very  satisfactory  results  for  three  reasons.  First,  the 
time  constant  of  the  system  includes  a  built-in  skew  in  the  data, 
especially  in  the  tails.  This  had  the  most  effect  upon  method  (1)  above. 
Second,  the  number  of  data  points  did  not  allow  an  accurate  peak  deter¬ 
mination.  Finally  (3),  the  overall  noise  level  caused  oscillations  in  the 
inverted  data.  This  is  very  evident  in  Fig.  20  which  shows  the  propa¬ 
gation  of  error  radially  inward. 


SECTION  VII 
SUMMARY 


Of  the  systems  utilized,  the  rotating  refractor  plate  synchronized 
to  the  oscilloscope  provided  the  most  flexible  system.  It  provided 
better  resolution  than  photographic  techniques  because  a  very  small 
exit  slit  could  be  used  (10  u).  Whenever  time  is  not  a  factor,  as  in  the 
case  of  continuously  running  facilities,  the  refractor  plate  technique 
can  provide  high  signal -to-noise  data. 

One  restriction  the  refractor  plate  imposed  was  dead  time.  Data 
were  only  taken  over  a  small  portion  of  the  plate  revolution,  0.  8  sec, 
and  the  total  time  between  scans  was  about  30  sec.  A  resetable  scanner 
might  provide  a  solution. 

The  photographic  techniques,  on  the  other  hand,  provided  a  large 
number  of  lines  from  which  to  work.  When  plates  are  to  be  exposed  for 
species  identification,  concentration,  or  other  such  spectroscopic 
studies,  very  little  extra  work  could  provide  velocities  as  well.  How¬ 
ever,  because  of  the  effects  of  emulsion  grain,  and  the  requirement  of 
spanning  relatively  large  distances  on  a  plate  to  reach  a  reference  line, 
this  method  is  less  precise  and  generally  more  time  consuming  than  the 
other  methods. 

Of  the  data  reduction  methods  described,  the  more  straightforward 
methods  of  measurement  and  peak  displacement  proved  to  be  the  most 
accurate  as  well  as  requiring  the  least  number  of  data  points.  The  line 
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Fig.  18  Abel  Inversion  with  Intensity  Monitor,  Input  Data 
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Station  or  Radius,  in. 


Fig.  20  Comparison  of  Input  and  Inverted  Data  at  Peak  Intensites 
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shape  methods  gave  fair  agreement  with  the  peak  displacement  methods 
but  usually  suffered  from  improper  choice  of  base  line  or  non-Gaussian 
effects. 

Application  of  an  Abel  inversion  should  follow  (1)  an  increase  in 
the  signal-to-noise  ratio  and  (2)  a  study  of  the  instrumental  broadening 
effects  on  the  observed  lines.  Methods  presently  exist  to  eliminate  in¬ 
strumental  effects  from  special  classes  of  line  shapes  (Ref.  6).  Much 
of  the  inaccuracy  associated  with  the  inversion  is  attributable  to  manual 
data  manipulation;  early  conversion  of  the  data  to  a  form  more  easily 
handled  by  a  digital  computer  would  benefit  the  system. 
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APPENDIX 

COMPUTERIZED  FORM  OF  ABEL  INVERSION 


G  LEVEL  LB  RA01NV  DATE  =  70355  15/34/53 

SUtJR*JUT  l«Mfc  RAD  I J J V I  ft »  C  ,  I-  )  • 

UlrtEftSIUi  tll)tF(l) 

0 ( MENS  l  Lift  A(25,25)  ,Aim,A2<2),A3i3),A4(4),A5(5),A6(6),  A7(7>, 

♦A618) , AVI 9), A10( 10>,A11 (11)  ,  A12C  12 ) , A 1 3 C 1 3 ) , A 1 4  1 14) ,A15I 15), 

*A16<  16),A17U  7)  ,A16C  18>,A19I 19) , A20I 20) , A2l (2 1 > , A22 C 22) ,A23(23), 
*A24C  24),A25(25) 

kwLl  VALtwCE  <  A<  If  l)fAll  l))f  IAU,2>,  A2ID)  ,(A|  1,3)  ,A3(1>  )t 
MAC  1,4),A4(  1>)  tCAC  1,5),A5C1)>  .tAll,6)  , Ab 1 1 > > , I  AC  1 ,  7),  A7I1)), 

*  I  At  1,8),  ABC  1) )  ,(A(  1,9) ,  A9C1))  »(  At  1«10>  v  A101 1)  )  f  (Ml  til) » AllllUi 

*  Utl,  l2),A12m>,tA(l,13),A13(l)),tAI  l,14),A14t  1)  ),(A(1,15), 

*A15( 1) i ,1 AC  1, 16) ,A16( 1) )« C AC  1*17) ,A17C1) ) ,1A( 1 ,18  >  »  Aldl 1)1* 

*t At  1,19), A 19 1  l))t(A(l»20)  fA20U))«IAIl,21)fA21Cl))f|A(l,22)f 
*A22( 1 )  )  *  (  A  ( 1 , 2 3  ) ,  A23{ 1) ) « I  AC  1 ,24 ) , A241 1) ) , I A( 1 ,25 ) , A25C1) ) 

DATA  A  1/  1  •  5  7080/ 

OATA  A2/ 2. 25565 ,2. 45674/ 

DATA  AV2-0605  1,2. 69599,3.09748/ 

DATA  A4/2.02691,2.2j719,3.1 0298 • 3*62649/ 

OATA  A 5/ 2. 01706*  2.  13023* 2.43399 *3*46 836 1 4 *08753/ 

DATA  A6/2.01 128, 2 •08348*2. 2 5572 ,2*62587,3*80085,4*50156/ 

DATA  A7/2.0U802 ,2* 05838,2* 17189 ,2*38535, 2*80879 *4. 10745 *4.68057/ 
DATA  A8/ 2 .00600, 2. 04323 ,2 *12445,2*267 17, 2 *51312,2 .98257 ,4.39322, 
*5.23219/ 

OATA  A9/2. 00466, 2. 033  34 ,2. 0946b  ,2 . 198  13. 2  *3  6.>90 ,2 .63750 , 3*  14795,  ' 
*4.66179, 5.50162/ 

OATA  A10/2. 003 72, 2. 02652, 2. 0745 7, 2. 15358, 2* 2 7468, 2*45988, 2. 75790, 

*  3.30579.4.91590,5.87259/ 

LiAT  A  A1 1/2.00304,2*02161,2*  Oo 03 5, 2*  12291 , 2.2 1602, 2* 35  19  2, 2.  55405 , 

*  2.87429, 3.45690,5.15763,6. 16790/ 

U AT  A  A 12/2. 0025 3, 2. 01 795,2. U49o9 ,2. 100 78 ,2 • 1 7499,2 . 2799 5, 2. 42875, 

*  2. 6463 1,2. 9868b ,3. 60200, 5.36861 ,6.44971/ 

DAT a  A13/2 .0J214, 2.01 515,2.04 195 ,2.06425 ,2. 14496,2*22899 , 2.34425,  ' 

*  2.50461,2. 7362^.3.09580,3. 7417 1 ,5. 61016 ,6. 71970/ 

DATA  A 14/2 .00 184, 2 .01 296, 2. 0357 9, 2. 07 153, 2 .12228, 2.  19 13 1, 2. 28 384,  ' 

*  2 .40830, 2 .57923,2. 82389, 3.2 0140, 3. 37*655, 5.82 332, 6.9  792  5/ 

DATA  A15/2 .00159, 2. 01122 ,2 .03090, 2.00154, 2 .10464,2* 16254, 2. 2 3 880 • 

*  2.33892, 2.47176,2.65246,2.90938, 3.30369, 4. 00699, 6 . 02899 , 7. 22950/ 
OATA  Alo/2 .00139, 2.00980,2. 02695,2.05352 ,2 .09 0o4, 2. 14000 ,2*20414,  “ 

*  2 .28682, 2. 39384, 2. 53441, 2. 7242 5, 2. 99277, 3 .40 349, 4. 13340 , 6. 22 790 , 

*  7.47136/ 

OATA  A1 7/2. 00 12 3, 2 .00 864, 2 .023 72 ,2 . 04699 ,2 .07932, 2 . 1 2 19 6,2. 1 7679, 

*  2. 2464/, 2. 33500, 2. 44836, 2. 59625, 2. 79479 , 3.0741 5, 3 .50043, 4. 2561 4,  ' 

*  6.42087,7.70564/ 

DA  I A  A18/2*  00109,2*  00  768,2. 02 10*«,2.  04  160, 2.07002,2. 1072  8 ,2*  1547  7, 

*  2. 2 l44o, 2 .2892 1,2 .38307, 2. 50234 ,2. 65662, 2. 06357, 3. 1536 5, 3* 59486, 

*  4.37550,6.60783,7.93300/ 

DATA  A19/2. 00098, 2 .00686, 2 .01 879, 2. 03710, 2. 06230, 2. 09516, 2. 13674, 

”  *  2 . 18855, 2 .25266,2 . 33 200, 2*43086 ,2. 55568,2.71667 ,  . . . 

*  2 .93115, 3.23137,3. 68698, 4. 49174, 6. 78985, 6. 15402/ 
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FUKTKAh  IV 
0024 


0025 


0026 


0027 


0028 


002V 


0030 

0031 

0032 

0033 

0034 

003b 

0036 

0037 

0036 

0039 

0040 


G  LEVEL  16  ftAOINV  DATE  ■  7035b  15/34/53 

Da  I A  A20/2.00U$8,2.006l7,2.0l669,2.03329,2.05980, 2.08502,2. 12178, 

*  *  2.16724,  2. 2229of  2 .201*1 ,2. 3  7473, 2  .*7826, 2.  *0833 , 2.77545, 2.99741, 

*  3.30740,3. 7769 3, 4 .60609, 6. 96 713,8.36 9*1/ 

DATA  A 2 1/2. 00079,2. 006b U, 2. 01 526, 2. 03005, 2. 05028, 2. 07644, 2. 10920, 

4  2.14946,2.19042,2.25773,2.32903,2.41727 , 2 . 52 520, 2 .6602 5, 2. 83324, 

*  3.06241, 3.38182,3.86 486, 4.71574, 7.14002, 8.57900/ 

DA 7*  A22/2. 000  7 2, 2. 00 507,2.01386, 2 .02726, 2 .04505, 2 .06962, 2. 09852, 

*  2.1 3 446, 2. 177U 8, 2. 23006, 2. 29269, 2 .36609,2 .45951, 2. 5 7163, 2. 7 1145, 

*  2.89006,3.12619,3.45473, 3.95069,4.02389,7.30882,8.78378/ 

DATA  m23/2. 00066, 2. 00463 ,2. 01 264, 2. 02484, 2. 04 146, 2. Oo 28 2, 2. 08936, 

*  2. 12 166, 2. 16048 ,2. 20682, 2. 26 197, 2. 32 7 7 1,2.4064 1,2 .50 141, 2. 61 752, 

*  2.76192, 2.94596,3. 18882,3.52621,4.03513,4.92967, 7.4738 3, 8. 98390/ 
DaTA  A24/ 2 .00060, 2.00424,2.01 158,2.02274 ,2.03790, 2 .05736,2.08145, 

*  2.1 1066,2.14560,2. 18 707, 2. 23609 »2.294U2 ,2. 36268, 2. 4445 1, 2. 54293, 

4  2. 66286, 2. 81166,3.00095, 3. 25033, 3.59634, 4. 11 768# 5. 03327, 7. 63527 , 

*  9.17465/ 

DATA  A 25/2 • 00056,2 • 00 390, 2.01 06 4, 2 .02089 ,2.03479, 2.05258, 2. 07456* 

*  2.10112,2.13275,2. 17012, 2.21404,2.26557 ,2.32612, 2. 39753, 2. 48234, 

*  2. 58403,2. 70764, 2. 86069, 3. 05506, >.310 78 ,3. 66518, 4. 19865, 5. 13479, 

*  7.79337,9.37132/ 

F  iS  *  r*  *  1 1 

Du  2  J=2 ,  N 

J  D  1=  JD  ♦  l 

M JO)=FN+F( JU) 

DU  3  K=JD1  ,N  .  .  '  ' 

3  CU6)  =  C(  J0J-A1  JD,K)*CIKi 

2  L  ( JOi  =C(  jU)/A(JD,JD)  ~  '  . .  ....  . - 

RETURN 

end 
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